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ABSTRACT 

This paper proposes a unified approach to joint adaptive parame- 
ter estimation and interference cancellation (IC) for direct sequence 
code-division-multiple-access (DS-CDMA) systems in multipath chan- 
nels. A unified framework is presented in which the IC problem is 
formulated as an optimization problem with extra degrees of free- 
dom of an IC parameter vector for each stage and user. We propose 
a joint optimization method for estimating the IC parameter vec- 
tor, the linear receiver filter front-end, and the channel along with 
minimum mean squared error (MMSE) expressions for the estima- 
tors. Based on the proposed joint optimization approach, we derive 
low-complexity stochastic gradient (SG) algorithms for estimating 
the desired parameters. Simulation results for the uplink of a syn- 
chronous DS-CDMA system show that the proposed methods sig- 
nificantly outperform the best known IC receivers. 

1. INTRODUCTION 

High data rate applications for future wireless systems require an 
ever-increasing sophistication and performance of receivers. In multi- 
user systems such as DS-CDMA, novel signal processing techniques 
are of crucial importance to enhance the capacity and the perfor- 
mance. The field of interference mitigation techniques has become 
an important and vibrant field since the pioneering work of Verdu 
m . The optimal multiuser detector has been proposed by Verdu in 
f2l, however, prohibitive complexity makes its deployment infeasi- 
ble and motivated the development of several suboptimal schemes 
that are amenable to implementation: The linear |3| and decision 
feedback |4| receivers, the successive interference canceler (SIC) 
(Sl . and the parallel interference canceler (PIC) IS). These receivers 
require the estimation of various parameters in order to carry out 
interference suppression. 

In most practical scenarios, the parameter estimation of the re- 
ceiver has to be computed adaptively in order to track the time- 
varying multipath channel conditions. Several adaptive receivers 
have been reported in I?] H] |9] and have proven to be very valu- 
able techniques for interference mitigation. Specifically for uplink 
scenarios, SIC IS'TO'TTl and PIC I'SIJH] [III [H [B] [IF| receivers, 
which are relatively simple and perform interference cancellation (IC) 
by sequentially or iteratively removing multiple access interference 



(MAI), are known to provide significant gains over RAKE and linear 
detectors. The works on SIC and PIC detectors present severe lim- 
itations with respect to the amount of interference to be estimated 
and cancelled in dispersive and dynamic environments. This is be- 
cause SIC and PIC detectors rely on an amplitude estimate of the 
parameters to be cancelled which affect the whole IC procedure. 

This work proposes a unified approach to joint adaptive IC de- 
tectors for DS-CDMA systems in frequency selective channels. A 
novel framework in which the IC problem is formulated as an op- 
timization problem of an IC parameter vector for each user and a 
unification of existing IC under the same model is described. A joint 
optimization method for estimating the IC parameter vector, the lin- 
ear receiver front-end filter and the channel parameters is proposed 
along with MMSE expressions. Based on the proposed joint opti- 
mization approach, we also present low complexity SG algorithms 
for estimating the desired parameters. Simulations for the uplink of 
a DS-CDMA system show that the proposed methods significantly 
outperform the best known IC receivers. 

The rest of this article is organized as follows: Section |2] de- 
scribes a synchronous DS-CDMA system model. In Section |3] we 
present the novel unified framework for IC. Section|4]details the pro- 
posed joint optimization approach and describes the MMSE expres- 
sions for channel estimation, receiver filter and IC parameter vector 
estimation. Section[5]is devoted to the proposed SG adaptive algo- 
rithms for joint IC and parameter estimation. Section |6]is dedicated 
to the simulations and to the discussion of the results, whereas the 
conclusions are drawn in Section [7] 

2. DS-CDMA SYSTEM MODEL 

Let us consider the uplink of a symbol synchronous quadrature phase- 
shift keying (QPSK) DS-CDMA system with K users, N chips per 
symbol and Lp propagation paths. Note that a synchronous model 
is assumed for simplicity, although it captures most of the features 
of asynchronous models with small to moderate delay spreads. The 
baseband signal transmitted by the fc-th active user to the base station 
is given by 



Xk{t) = Afc ^ hk[i]sk{t~iT), 



(1) 
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where hk[i\ £ {±1 ± j} with j^ — —1 denotes the i-th symbol 
for user k, the real-valued spreading waveform and the amplitude 
associated with user k are Sk(t) and Ak, respectively. The spreading 
waveforms are expressed by s^ (t) = X]i=i '^k [i\4>{t ~ iTc), where 



ak[i] G {lizl/VN}, (j)(t) is the chip waveform, Tc is the chip du- 
ration and A^ = T/Tc is the processing gain. Assuming that the 
receiver is synchronized with the main path, the complex envelope 
of the coherently demodulated composite received signal is 



K Lp 



'"W ^ 1] 5Z hk,i{t)xk{t - rk,i) + n{t), 



(2) 



where hk,i{t) and Tk,i are, respectively, the channel coefficient and 
the delay associated with the Z-th path and the fc-th user. Assuming 
that Tk.i — ITc, the channel is constant during each packet transmis- 
sion, Lp < N and the spreading codes are repeated from symbol 
to symbol, the received signal r{t) after filtering by a chip-pulse 
matched filter and sampled at chip rate yields the Af -dimensional 
received vector 



r[i\ = ^ Afc (bk[i - i]c5; + bk[i\Ck 



(3) 



+ bk[i + l]Cl)hk + n[i 



where M = N + Lp — 1, n[i] = [iT-ili] ■ ■ ■ nM[i]]'^ is the com- 
plex Gaussian noise vector with E[n[i]n^[i]] — a^I, where (■)"^ 
and {■)^ denote transpose and Hermitian transpose, respectively, 
E[-] stands for expected value, the amplitude of user k is Ak, s^ = 
[afe(l) . . . ak{N)]^ is the signature sequence for the fc-th user, the 
AI X Lp constraint matrices C^, Ck, Cf, that contains one-chip 
shifted versions of the signature sequence for user k and the Lp x 1 
vector hfc with the discrete-time multipath components are described 
by 
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(5) 

where the matrices C^ and Cf account for the intersymbol interfer- 
ence from the previous and subsequent symbols, respectively. The 
subscript [m : q,j : p] denotes the range of elements of a given 
matrix used. 

3. UNIFIED FRAMEWORK FOR INTERFERENCE 
CANCELLATION 

In this section, we present a unified framework for IC based on the 
formulation of the problem as the optimization of an IC parameter 
vector for each user. Prior works on IC are heavily based on the 
estimation of an amplitude of a user (or an amount of interference) to 
be reconstructed and cancelled. This approach is very sensitive and 
tends to be inaccurate in multipath scenarios. Unlike prior works, 
the proposed approach trades off an amplitude estimate against a 
parameter vector estimate for each stage, incorporating more degrees 
of freedom to the problem and which reveals to be significantly more 
effective than the best known methods. With this new formulation, 
we provide a unifying treatment of IC schemes in what follows. 



Let us first consider a conventional IC approach, where the cen- 
tralized processor aims to reconstruct the detected data and subtract 
them from the M x 1 received data 



Ft t = r 



M-^i7(b7[i^i]c^ + fe™wc, 



where F7[i 



jeg 
+ b™[i + l]C,= )h™[i] 

jeg 



= ibf[i 



(6) 



M X Lp matrix with the signature code and symbol estimates of user 
j at the IC stage m, Q — {Qi, Q2, ■ ■ • , Gp} denotes the group of 
users with P entries to be reconstructed and subtracted, whereas the 
subscript denotes the fcth user. In the conventional approach outlined 
in l[6}, the goal is to compute the estimates A™ , F™ [i] and h™ [i] and 
this often leads to inaccurate IC for systems in dispersive channels. 
The novel approach we detail here corresponds to a mathemat- 
ical reformulation of (|6) and the introduction of an IC parameter 
vector A™ = [AJI^^ A™2 ■ • ■ ^T,p] for each stage m and user k, 
as expressed by 



jeg 



■ tF™ 



\i]h" 



(7) 
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where the MxP matrix with signature codes, symbols and channels 
estimates of the group of users Q is given by 



m-r-tm. 



Dg[i\=C^H"'B 

Ct = Cgj Cg, 



(8) 
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(10) 
(11) 



B^H = diag^^fog'Jj], bgjj], . . . , Ug^ 

The matrix Dg'[i] corresponds to the reconstructed data of users 
which belong to group Q and it is a function of the M x [P ■ Lp) 
code matrices C^, Ct, Ct, the P x P user's symbol matrix B^fi] 
and the {P ■ Lp) x P channel matrix 'H"^ with a block diagonal 
structure, with all the parameters of the users to be reconstructed 
and cancelled. 

The mathematical framework detailed in (|7j can be actually used 
to describe IC schemes which perform SIC and PIC as particular 
cases. For instance, if the designer chooses to detect the users ac- 
cording to a decreasing power ordering, we obtain the following SIC 
approach 

fc-i 
rT[i] - r[i\ - ^ A™fcF™[i]h™[i] 

= r[i] -Dg^_Jj]A™[i], m = l 

where Gk-i ~ {Gi, G2, ■ ■ ■ , Qk-i} denotes the group of users to 
be reconstructed according to the SIC approach (decreasing power 
order). Note that there is only one stage (m — 1) for SIC, user 
k is detected and the previously detected users are regenerated and 
subtracted from the received data r[i] and this is repeated for the 
remaining users. 
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Fig. 1. Proposed receiver structure. 



Another detection strategy whiicli can be carried out as a partic- 
ular case of tlie framework in (|7) is tlie multistage detection or PIC. 
Following the proposed detection scheme users are detected on the 
basis of the following structure 



Tfc W 
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(13) 



Dm r '1 \ '^ 



m = 1, 2, 



where k represents the desired user to be detected, m is the stage, 
and Qk-1 contains all but the desired fcth user. 

4. JOINT INTERFERENCE CANCELLATION AND 
PARAMETER ESTIMATION METHOD 

In this section, we present a novel strategy for joint interference can- 
cellation and parameters estimation based on the formulation given 
in the previous section. The idea is to consider the problems of re- 
ceiver filter, channel, and IC parameter vector estimation jointly and 
devise MMSE expressions to solve it. Specifically, we consider here 
that a linear receiver filter is employed at the front-end of the IC de- 
tector. The proposed receiver structure is showed in Figure [T] Let us 
first consider the following cost functions 



Ji(wrM)=E[l6fe[j 



Wfe Wrft [i]\ 



(14) 



:[\\FT\^]hT\^]-v[l]+-D^\^]\Tm']■ 

(15) 

By minimizing ( I14t with respect to the receiver linear filter w™ [i] 
we obtain the Wiener-Hopf-like expressions 



Rri[j]p6. 



(16) 



where Rr^ [i 



E ri: 



is the M X M covariance matrix 



and Pb;, [i] = E[6^[i]r™[i]] is the M x 1 cross-correlation vector. 
By taking the gradient terms of ( I15t with respect to the IC parameter 
vector A^" [i] and equating them to zero we get 



-^rw 



Rd^Hpf^ 



(17) 



where Rd^H = E[Dg '^[j]D5'[i]] is the PxP covariance matrix 

and pf" [i] = E[D™"^H(r[i] - F^TMhlT W)] is the P x 1 cross- 
correlation vector. By minimizing (115) with regard to the channel 



estimate h™ [i] we obtain the last system of linear equations 



fe 



r; 



mpo" 



(18) 



where RF"«[i] = E[F™'^[i]F™[i]] is the Lp x Lp covariance ma- 
trix and PDg'M = E[F™-^[i](r[i] - Dg'[i]A^[i])] is the P x 1 
cross-correlation vector The expressions in (|16t - (|181 are not closed- 
form ones as the interference cancellation and the receiver linear fil- 
ter parameters depend on the channel and vice-versa. This means 
that l|16t-l|18t have to be iterated in order to seek a solution to the 
optimization problem. The detected symbols are obtained as fol- 
lows: 

k[i] = sgn[m{xk\i])] + j[Q{xk[i])], (19) 

where Xk[i] ~ wj."' [j]r™[j] and sgn[ • ] is the signum function. 
The expressions in JI6t -(II8t require a computational complexity of 
0{M^), 0(P3) and O(L^) for the estimation of w^[i], A^[i] and 
h™[i], respectively. This cubic cost is due to the required matrix 
inversions. Note that the optimization of J2(A™[i], h^[i]) in l |15t 
seeks to eliminate MAI, whereas the optimization of Ji (w™ [i] ) in 
l ll5t aims to suppress the residual MAI and intersymbol interference. 
However, a problem is that the statistics are not known a priori and 
the designer has to estimate them. In what follows, we seek adap- 
tive SG solutions based on the joint optimization of the parameters 
formulated in ( II4t and dlSt for estimating these parameters and al- 
lowing amenable implementation. 

5. ADAPTIVE ESTIMATION ALGORITHMS 

In this section, we propose adaptive estimation algorithms based on 
the joint optimization problems stated in ( I14t and ( |I5t . In order 
to develop SG adaptive algorithms, we compute the instantaneous 
gradients of ( I14t with respect to w™ [i], and of J15t with respect to 
A™[i] andh^H as follows 
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: Fr'^W(FrWhrM - rW + D^WArH) (22) 
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■ k l«Jefc [i\ 

By using the well-known gradient-descent rules and introducing step 
sizes, we devise the following SG adaptive estimators: 

w^ [J + 1] = w^ [i] + Mii' efe "" [i] r™ [i] , (23) 

XT[i + 1] = ArW - MAD^'^HerW, (24) 

K [i + l]^ K W 't^nK'" Wer [i\ , (25) 



where /^u,, fix, and Hh are the step sizes for the adaptive SG estima- 
tors for the receiver filter, the IC parameter vector, and the channel, 
respectively, and scalar and vector errors are given by 



N=16, L = 9, K=8 users, Average ^/N|j=15 c 



J r -1 m,H T •! n 



(26) 



e^M = FrWhrW - rW + D^HArW (27) 

The adaptive algorithms described above should be iterated in order 
to converge to a solution since the parameters estimated by them 
are inter-dependent. The complexity of the algorithms in (I23t -ll27t 
is 0{M), 0{PM) and 0{MLp) for the estimation of the receiver 
filter w™ [i], the IC parameter vector A™ [i], and the channel h^ [i], 
respectively. 

6. SIMULATIONS 

In this section, we evaluate the bit error rate (BER) performance 
of the proposed joint interference cancellation and parameter esti- 
mation algorithms. We compare the proposed algorithms with the 
best known methods of interference mitigation, namely, the linear, 
the SIC, and the PIC detectors using SG algorithms as the estima- 
tion procedure. The DS-CDMA system employs randomly gener- 
ated spreading sequences of length A'^ = 16. The channels h^ = 
[hk,o hk,i ■ ■ ■ hk^Lp-i]"^ are modeled by a tapped-delay-line with 
Lp = 9 taps and are normalized such that X]i=o 1 ^; I ^ = 1- Specif- 
ically, of the Lp = 9 taps, there are only 3 non-zero paths with 
complex random gains, whose real and imaginary parts are gener- 
ated in each transmitted packet by uniform continuous random vari- 
ables (r.v.) between —1 and 1. The spacing between the paths in 
each packet is given by a discrete uniform r.v. between 1 and 3 chips. 
The system has a power distribution among the users for each trial 
that follows a log-normal distribution with associated standard devi- 
ation of 3 dB, the performance is shown in terms of average Et/No, 
and all curves are averaged over 100 runs. The packets used in the 
simulations have 1500 symbols and the training sequences have 150 
symbols. After the training sequences, the receivers are switched 
to decision-directed mode. The receiver filters used in the simula- 
tions have M — N + Lp — 1 = 24 taps. The step sizes of all 
algorithms are optimized for each situation to ensure the best BER 
performance for each packet. All the receivers considered in this 
part employ a linear receiver as the front-end. When there is no IC 
after the linear front-end, the structure corresponds to a conventional 
adaptive linear receiver f7l-f9l. The proposed jointly optimized (JO) 
algorithms with SIC and PIC receivers are denoted JO-SIC and JO- 
PIC, respectively. For the PIC receivers we used m = 3 stages and 
the amplitude estimation based on the estimate at the output of lin- 
ear receiver front-end [16 1 and for the proposed JO-PIC we also use 
m — 3 stages. For the existing SIC receivers llTolllII we used the 
amplitude estimation algorithm reported in Bill . 

In the first experiment, shown in Figure |2] we assess the BER 
convergence performance of the proposed adaptive estimation algo- 
rithms and receiver structures. The results indicate that the proposed 
JO-SIC and JO-PIC receivers have the best performance among the 
compared structures. The JO-SIC slightly outperforms the JO-PIC, 
which is followed by the SIC, the PIC, and the linear receivers. In 
particular, we notice that the convergence performance of the JO-SIC 
and JO-PIC is significantly superior to the remaining approaches. 
This is because the IC is substantially more accurate than the IC car- 
ried out by the existing SIC and PIC approaches. 

The second experiment, depicted in Figure [5] evaluates the MSE 
performance of the channel estimators. The channel estimators of 
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Fig. 2. BER performance versus number of received symbols. 



IC schemes clearly benefit by the cancellation process as compared 
to the linear estimator without IC. The use of IC for improving the 
performance is particularly relevant for the proposed JO-SIC and 
JO-PIC, which achieve the best performance. 

The BER versus Et/No and number of users is illustrated in 
Figure |4] The curves indicate that the best performance is obtained 
by the proposed JO-SIC and JO-PIC receiver and algorithms. The 
plots show that the proposed JO-SIC and JO-PIC receivers can save 
up to 4 dB in Eb/No for the same BER as compared with the lin- 
ear receiver, and up to 2.5 dB as compared with existing SIC and 
PIC detectors. In terms of system capacity, the proposed detectors 
and algorithms provide a substantial capacity improvement and can 
accommodate up to 50% more users for this small system. 

7. CONCLUSIONS 

This work proposed a unified framework for IC in DS-CDMA sys- 
tems, which formulates the IC problem as the optimization of an IC 
parameter vector. A joint optimization method for estimating the IC 
parameter vector, the receiver linear front-end filter and the chan- 
nel parameters along with MMSE expressions was also presented. 
Low-complexity adaptive estimation algorithms were developed for 
jointly estimating the desired parameters. The results for the uplink 
of a synchronous DS-CDMA system show that the proposed meth- 
ods significantly outperform the best known IC receivers. 
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